Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 19 May 2009 (MN WI^ stylo file v2.2) 



0^ 
O 
O 



O 

U 

6 



> 

cn 

in 
o 
a^ 
o 



Constrained simulations of the local universe: II. The 
nature of the local Hubble flow 

Luis A. Martinez- Vaquero^, Gustavo Yepes^, Yehuda Hoffman^, Stefan Gottlober^, 
Mira Sivan^ 

^ Grupo de Astrofisica, Universidad Autonoma de Madrid, Madrid E-280049, Spain 

^ Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 

^ Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14-4-82 Potsdam, Germany 

19 May 2009 

ABSTRACT 

Using a suite of N-body simulations in different Cold Dark Matter (CDM) sce- 
narios, with cosmological constant (ACDM) and without (OCDM, SCDM), we study 
the Hubble flow {(Th) in Local Volumes (LV) around Local Group (LG) like objects 
found in these simulations, and compare the numerical results with the most recent 
observations. We show that ACDM and OCDM models exhibit the same behavior of 
an- Hence, we demonstrate that the observed coldness of the Hubble flow is not likely 
to be a manifestation of the dark energy, contrary to previous claims. The coldness 
does not constitute a problem by itself but it poses a problem to the standard ACDM 
model only if the mean density within the Local Volume is greater than twice the 
mean matter cosmic density. The lack of blueshifted galaxies in the LV, outside of the 
LG can be considered as another manifestation of the coldness of the flow. Finally, we 
show that the main dynamical parameter that affects the coldness of the flow is the 
relative isolation of the LG, and the absence of nearby Milky Way like objects within 
a distance of about 3 Mpc. 
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1 INTRODUCTION 

The neighborhood of the Local Group (LG) is often de- 
scribed as being cold. This attribute implies that the dis- 
persion of the radial velocities of galaxies from a pure Hub- 
ble flow is s mall, and the ' s malln ess' amounts to less than 
100 km s~^. ISandage et al.l (|l972l ) studied local departures 
from a uniform Hubble flow and could only pu t upper 
limits on such departures. This led ISandage et al.l to con- 
clude tha t (?o S:^ (where go is t h e dece leration parameter). 
Later on ISandage fc TammannI (|l975h estimated that the 
upper limit to the mean random motion of field galaxies is 
< 50 km s~^. These early findings of Sandage, Tammann 
and their collaborators have be en corroborated and vigor- 
ously improved by many others. iKarachentsev et al.l (|2003l ) 
estimated the radial peculiar velocity dispersion of all galax- 
ies within 5.5 Mpc to be 85 km s~^. This value drops down 
to 41km s~^ if members of galaxy groups are removed and 
distance erros are t aken into account. U s ing a newer set of 
Karachentsev's data lTikhonov fc KlvpinI (|2008l ) found a ve- 
locity dispersion of 97km s~^, within 7 Mpc, which reduces 
after correction for apex motion and d istance errors only 
slightly to 84 km s~^. lMaccio et aEI (|2005l ) compiled the data 



from three different sources: The Cephei d-based distance 
meas urements of the HST Key Project ijFreedman et al.l 
l200ll ). distance estim ates based on the surface brightness 
(|Tonrv et all l200ll ) and TuUy-Fisher 



fluctuations method 



distances (|Tu11v et al1ll992l ). They fitted the data by an 



1 ± 20 km s"^ X (i?/7Mpc), where an is a measure of the 
dispersion of the radial velocities around a pure Hubble flow 
of galaxies within a sphere of radius R (a thorough discus- 
sion of the various estimates of gh is given below) . 

The observational evidences for a local cold Hubble flow 
seem to be indisputable. Yet, the question arises as to why 
a an of the order of a few tens of km s~^ is labeled as 
'cold'. Namely, by what standard is it cold?. Rich clusters of 
galaxies provide the first and the most robust evidence for a 
departure from a pure Hubble flow, with a dispersion of pe- 
culiar velocities of up to ~ 10^ km s~^. Compared with the 
rich clusters the neighborhood of the LG is definitely cold. A 
statistical estimate is also given by the pair weighted velocity 
dispersion {(J12) which was measured from the CfA redshift 
survey to be a\2(r = lh~^Mpc) = 340 ± 40km s~^ (Davis 
& Peebles 1983). Another more robust measure of the devi- 
ation from the Hubble flow is provided by the g\ statistics 
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which measures the one-dimensional rms peculiar-velocity 
dispersion of galaxies relative to their neighbors within a 
projected radius of 2h~^Mpc (Davis, Miller & White 1997). 
These authors found oi = 95 ± 16 km s"'^ (for the IRAS 
survey) and 130 ±15 km s~^ for the UGC catalog. The mea- 
sured (Ti is indeed mu ch 'hotter' than the an = 25 km s~^ 
within R = 3.0 Mpc (|Karachentsev et all [2OO9I ). So, with 
regard to the other measure of the dispersion of peculiar 
velocities the immediate neighborhood of the LG is indeed 
very cold. However, one should recall that the cti and G12 
measures consider all galaxies in a given survey. The cth 
considered here, on the other hand, refers to one particular 
object, namely the LG, that resides in a dynamically quiet 
environment. It should be emphasized that the LG is not an 
atypical object, yet the LG environment is not a representa- 
tive one for Milky Way like objects. This implies that, when 
addressing the issue of 'coldness' by means of simulations, 
the selection of the LG-like objects should be carefully done. 

The standard model of cosmology consists of a flat 
Friedmann universe whose mass energy is dominated by a 
cosmological constant (A) and cold dark matter (labeled as 
ACDM). It has been recently stated that the cosmologi- 
cal constant, or its generalization the Dark Energy (DE), 
manifests itself in the dynamics of the very loc al universe 
jBarvshev et aLlbOOllTchernin et alll2004[2007h . These au- 
thors argued that the coldness of the local flow is a mani- 
festation of t he existence o f the d ark energy. This has been 
supported bv lMaccio et all (|2005l ') who analysed a set of N- 
body simulations and concluded that indeed ".. [their] re- 
sults provide new, independent evidence for the presence of 
dark energy on scales of a few megaparsecs" . These results, 
if correct, would have provided an independent corrobora- 
tion to the DE component whose existence is otherwise in- 
ferred from observations of distant objects and on very large 
scales of the Universe. T he claims of^Ch crmn et al.l (j2007l ) 
have been challenged by iHoffman et al.l (|2008l 'l. who ana- 
lyzed a suit of constrained and unconstrained CDM simu- 
lations, identified LG-like objec ts and studied the flo w field 
around these. The main result oflHoffman ct al.' ('2OO8I) inval- 
idates the model advocated by[ciicrnin c t al. ( 200 7) ) and its 
basic predictions. The claims made bv iMaccio et all (|2005l l 
have remained unchallenged till present. 

There are two main issues this paper addresses. The 
claim that the coldness of the local flow is induced by the 
dark energy, and the fact that it is strongly contradicted by 
iHoffman et al.l (|2008l l. motivated us to revisit the problem. 
Then, there is the general problem of the nature of the local 
flow. In particular the question of how the various character- 
istics of the local neighborhood are afl'ecting the local flow 
is addressed here. 

The present p aper adopts the methodology used b y 
iMaccio et aP lj2005l ) and later on bv IHoffman et all (|2008h . 
Namely, a set of numerical simulations is performed, ensem- 
bles of LG-like objects are constructed and oh is calculated 
around these mock objects in the same way as the Hub- 
ble flow around the actual LG is analyzed. This is further 
extended here. The cosmological models studied here are 
the presently accepted model of cosmology with a cosmo- 
logical constant, ACDM, the open CDM (OCDM) and the 
outdated Einstein-De Siter SCDM models (where S stands 
for the standard of the years of the 1990s and CDM stands 
for cold dark matter). Most of the simulations studied here 



are constrained and thereby are constructed to reproduce 
within the simulation box the observed Large Scale Struc- 
ture (LSS) on scales larger than ~ 5h~^Mpc. A thorough 
study of the nature of the local Hubble flow is done by study- 
ing the dependence of an on a multitude of factors, including 
the cosmological parameters, the mass resolution, the very 
local environment and the structure on larger scales. The 
control over the very small scales is achieved by varying the 
selection rules for the LG-like objects and the control over 
the large scale environment is done by performing the con- 
strained simulations. A somewhat different selection of mock 
objects is obtained by constraining them to have no nearby 
blueshifted galaxies. 

The structure of the paper is as follows. A summary 
of the CNG data and examination of the different ways of 
calculating cth are given in fJS] §3] describes the N-body sim- 
ulations analyzed here and a description of the selection of 
the LG-like objects is presented i n j}4l The analysis of oh is 
conducted first along the lines of lMacci6 et all (|2005D (©. 
Then the dependence of oh on the various characteristics 
of the LG is shown in §S] and the possible role of mass res- 
olution, in SJT] The lack of nearby blueshifted galaxies as a 
manifestation of the coldness of the local flow is analyzed 
in ^ The paper concludes with a summary and a general 
discussion (Sj9|. 



2 OBSERVATIONS 

The Catalog of Neighboring Galaxies (jKarachentsev et al.l 
|2004| . CNG) provides the latest and most comprehensive sur- 
vey of the velocities of nearby galaxies. The catalog includes 
the distances and radial velocities of more than 400 galaxies, 
roughly 300 of which are in the Local Volume (LV), defined 
by a s phere of i? = 7 Mpc about the LG. lTikhonov fc KlvpinI 
(|2008l ) have recently compiled the catalog and calculated an 
of all the CNG galaxies within t he LV. Distances are typ- 
ically measured with 10% errors. iTikhonov fc KlvpinI stud- 
ied the recent, yet unpublished version of the catalog. That 
compilation of the data is taken here as the observational 
reference value, against which our numerical results are to 
be compared. 

The l ocal flow is studied here within the framework 
deflned bv ITikhonov fc KlvpinI (|2008l ) . Their study ignores 
non-linear motions within the LG and therefore all galax- 
ies closer than 1.0 Mpc are excluded from the analysis. The 
flow is analyzed by means of measuring the dispersion of the 
radial velocities about the Hubble flow of all galaxies found 
within the LV, which is redefined to correspond to the shell 
of (1 — 7) Mpc. The choice to consider all galaxies regard- 
less of their possible membership in galaxy groups is driven 
by the principle of keeping the selection to be as clear and 
simple as possible. This avoids the issue of defining a galaxy 
group catalog and simplifies the comparison of observations 
and simulations. 

There are different ways of defining the scatter of the 
line-of-sight (LOS) velocities from a pure Hubble, oh- Con- 
sider a catalog of galaxies whose dist ances and radial ve- 
locitie s are measured, {(r^, Di)}i=i,...,jv. ITikhonov fc KlvpinI 
(|2008l ) estimated the apex motion of the observer, namely 
the LG, with respect to the sample of galaxies in the cata- 
log and subtracted the apex velocity from the observed LOS 
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velocities. The dispersion an is then calculated as the root 
mean square (RMS) with respect to the Hubble flow, where 
the global va l ue of the Hubble constant (Ho) is assumed. 
iMaccio et al.l (|2005l ) defined an by the standard deviation 
(SDV) of the residual LOS velocities. This is equivalent to 
measuring the RMS of the residual velocities from the lo- 
cal Hubble flow, namely using the locally determined Hub- 
ble constant from the data itself (Hioc)- Yet, a simpler ap- 
proach is to associate an with the RMS of the residual from 
the global Hubble flow. This is the approach adopted here, 
namely an is calculated by: 
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and the sum extends over the A*' galaxies within the appro- 
priate distance cut of the catalog. Hereafter these methods 
are referred to as the APEX, SDV and RMS. 

One can argue for the merits and disadvantages of the 
different methods and they all can be used as long as obser- 
vations and simulations are analyzed consistently. Still, we 
argue here that the use of the dispersion around the global 
Hubble flow is the preferred way. The motivation for mea- 
suring an stems from the information it provides on the 
deviations from a pure Hubble flow. By using the locally 
determined Hubble constant the 'breathing' (i.e. isotropic) 
mode of the flow is absorbed in Hioc, and so the resulting 
an underestimates the actual deviation from a pure Hubble 
flow. In the case of the apex motion the velocity of the LG, 
that hosts the observer, is treated as if it is an external ef- 
fect that is not related to the perturbed flow one is trying 
to study. It is solved in a manner which does not depend 
on the nature of the perturbation field. The inclusion of 
the apex motion of the observers provides a better measure 
of the deviation from a pure Hubble flow. Yet, it is shown 
here that the differences between the different measures are 
small. The different measures of an have been applied to 
an ensemble of LG-like objects drawn from one of the high 
resolution simulations studied here (ACDMhr, see ^ for a 
detailed description). Figure [T] presents the results obtained 
for the different estimators of an - The main result is that the 
differences between the different methods are much smaller 
than the scatter around the median of calculated values of 
an- 



Table 2 of lTikhonov fc KlvpinI lj2008l ) provides the RMS 
and APEX estimates of an for the CNG galaxies for various 
radial cuts. These estimates are corrected for distance errors. 
The representative number used here is the RMS estimation 
for all galaxies within (1 — 7) Mpc of an = 90.4km s~^. 



3 SIMULATIONS 

A suit of constrained and unconstrained ACDM (flat, A and 
CDM), OCDM (open, CDM only) and SCDM (flat, CDM 
only) low and high resolution simulations have been per- 
formed. Table [T] lists the names and the parameters of the 
simulations. The five low resolution simulations (A'^ = 256^) 
are th e sa me ones descr i bed i n iMartinez-Vaauero et al.l 
|2003) and iHoffman et all (|2008l '). A computational box of 
L = 64h~^Mpc (where h is Hubble's constant in units of 
100km s"'^ Mpc~^) and the WMAPl cosmolo gical parame- 
ters are assumed for the ACDM simulations (jSpergel et al.l 
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Figure 1. The median of au (lines with points) and its standard 
deviation (thin Unes) for the ACDMhr simulation taking all the 
DM halos (solid lines) and only those with a circular velocity 
higher than 35km (dotted lines) in the distance range of (1 — 
_R) Mpc computed using the APEX (triangles), SDV (squares) 
and RMS (circles) methods. 



l2003h . The OCDM simulations are identical to the ACDM 
ones but without the A term, hence they correspond to 
an open universe. The SCDM is used here just as a ex- 
treme model for which the Hubble flow is expected to be 
much hotter than all other models, despite the fact that is 
a cosmological model which stands in marked disagreement 
with a wide range of cosmological observations. The five 
low resolution, constrained and unconstrained, simulations 
are all based on the same realization of the Gaussian ran- 
dom field. The ACDMhr and ACDM160 are high resolution 
(A'' — 1024^) simulations in computational boxes of L = 64 
and 160h~^Mpc respectively. Both are performed with the 
WMAP3 cosmological parameters (|Spergel et al.|[2"007^ . 

Five of the simulations are constrained ones, namely the 
initial conditions of these are set by constrained realizations 
of Gaussian fields. Observational data of the nearby universe 
is used as constraints on the initial conditions and thereby 
the resulting simulations reproduce the observed LSS. The 
basic algorithm of constraining Gaussian random fields is the 
one developed by Hoffman & Ribak (1991). The implemen- 
tation of the algorithm to observational data and a descrip- 
tion of the constru ction of co nstrained simula tions was de- 
scribe d at length in lKravtsov ct al. ( 20^ ') and lKlvpin et all 
(j2003f l. A brief description of the constraining data is given 
here. Two different observational data is used to set up the 
initial conditions. The first is made o f radial velocities o f 
galax ies drawn from th e MAR K HI j WilUck ct al.' 'l997|), 
SBF (jTonrv et al.ll200lf ) and the lKarachcntscv ct al._ (20M) 
catalogs. Peculiar velocities are less affected by non-linear 
effects and are used as constr aints as if they were linear 
quantities (jZaroubi et al.lE999[ l. The other constraints are 
obtained from the catalog of nearby X-ray selected clusters 
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of galaxies (|Reiprich fc Bohringeij|2002l ). Given the virial pa- 
rameters of a cluster and assuming the spherical top-hat 
model one can derive the linear overdensity of the cluster. 
The estimated linear overdensity is imposed on the mass 
scale of the cluster as a constraint. For the CDM cosmogo- 
nies the data used here c onstrains the simul ations on scales 
larger than fa 5 h"^Mpc (|Klvpin et al.ll2003l ). It follows that 
the main features that characterize the local universe, such 
as the Local Supercluster, Virgo cluster, Coma cluster and 
Great attractor (in the large box), are all reproduced by the 
simulations. The small scale structure is hardly affected by 
the constraints and is essentially random. 

We hav e used the pa rallel TREEPM N-body code 
GADGET2 (|Springell [ioosi ) to run these simulations. An 
uniform mesh of 512^ grid points was used to compute the 
long-range gravitational force by means of the Particle- Mesh 
algorithm. A constant comoving Plummer equivalent grav- 
itational smoothing scale of 20h~^kpc was set at high red- 
shift and we changed it to 5h~^kpc kpc physical scale since 
z=3 till z=0. The simulations were started at z=50 in all 
cases. We employed a variety of parallel computer architec- 
tures (SGI-ALTIX, IBM-SP4, Opteron-clusters) during the 
course of this work. Using 16 processors simultaneously, we 
completed one run in about 2 cpu days. 

We analyzed also two simulations with higher mass 
resolution (1024^ particles). The simulation with box size 
64h"^Mpc has been also evolved using GADGET2. A grid 
of 1024"^ was used in this case for the computation of the 
PM forces. Gravitational smoothing was set to 1.6h~^kpc 
at high redshift and 0.8 h~^kpc at lower redshift. The ini- 
tial conditions were set up at z = 100 for this run. The 
simulation of box size 160h"^Mpc has been calculated by 
the MPI version of the Adaptive Refineme nt Tree (ART) 
code described in lGottlober fc KlvpinI (|2008l ). Initial condi- 
tions were set up at z = 30 and the maximum number of 
refinement levels were set to 9, which translates into a maxi- 
mum spatial resolution of 1.2 h~^kpc. T hese two high resolu- 
tion s imulations were also analyzed bv lTikhonov fc KlvpinI 
(|2008D . 

We conclude this section with a general remark. The 
parameters of the simulations are defined by using the h^^ 
scaling. Yet, the comparison of the results with observations, 
and in particular the distance cuts, is made after an Hq = 
73 km s~^ Mpc~^ is assumed and distances are expressed in 
units of "real" Mpc. 



4 SELECTION OF LG CANDIDATES 

DM halos were found in simulations using two object find- 
ing methods: The B ound Density Maxima (BDM) algorithm 
(|Klvpin et al.lll999l l is based on finding local center of mass 
in spheres of variable radius starting from randomly se- 
lected particles i n the simulation. The AMIGA Halo Finder 
iGill et al.ll2004l ). on the contrary, finds local density max- 
ima from an adaptive mesh hierachy. In both cases, an it- 
erative procedure to find local centre of mass from density 
maxima is used. Particles that are gravitationally unbound 
to the halo potential are also removed. Only halos with a 
mass higher than 2.6 x 10^° (ACDM, ACDMu, OCDM and 
OCDMu), 8.7 X 10^° (SCDM), 5.0 x 10* (ACDMhr) and 
4.0 X l(fMQ (ACDM160) are considered. For the work re- 



ported here we have used the halo catalogues obtained by the 
AMIGA code, except for the ACDM160 simulation where 
we have used the BDM halo catalogue. In any case, we have 
checked that the results are independent of the halo finder. 

One of the main a ims of this paper is to perform a de- 
tailed comparison with lMaccio et al" I (|2005h results, so as to 
be able to challenge their claim on the role of the DE in the 
local dynamics. In order t o meet this go al the selection cri- 
teria of LG-like objects of iMaccio et al.l have been followed 
to the letter. Going beyond the comparison with that paper 
the selection rules are reexamined and some of which are 
revised. Here, these criteria are expressed in physical units 
without the scaling. The LG objects obey the following 
selection: 

(i) The groups contain two MW and M31 like DM halos 
with maximum circular velocity in the range of 125 ^ Vc ^ 
270 km/s. 

(ii) The two major DM halos are separated by no more 
than 1 Mpc. 

(iii) The relative radial velocity of the two main halos is 
negative. 

(iv) There are no objects with maximum circular velocity 
higher than MW and M31 candidates within a distance of 
3 Mpc. 

(v) The group resides within a distance of 7 to 17 Mpc 
from one and only one Virgo like halo of 500 ^ 14 ^ 
1500 km/s. No Virgo like halos can appear within a distance 
of 7 Mpc. 

The selected LG-like objects obeying all the above rules 
are defined as Pairs, for the fact that they are dominated 
by two MW and M31 like DM halos. 

A careful examination of the dynamics of LG-like ob- 
jects suggests the dynamics of the Hubble flow might not 
depend strongly on whether they are dominated by two, 
almost e qual, massive object s or by rather one massive ob - 
jects (see lMaccio et al.ll2005l and lTikhonov fc Klvpinll2008l '). 
To test this idea we construct the Singles ensemble of LG- 
like objects which obey all the criteria presented above apart 
from the first one. Namely, for the Singles objects we look 
for individual halos, whose mass are similar to MW and M31 
together, namely 10^^ M^ir ^ 2 x 10^^ h"^MQ. The LV 
around such objects is studied. 

In Table [2j the number of LG candidates found follow- 
ing both criteria is shown. Higher number of LG-like objects 
are obtained using Singles criterion becuse it is less restric- 
tive than the Pairs one. One should note that there is not a 
one-to-one correspondence of mock LGs in the ACDM and 
OCDM simulations, even if they have the same computa- 
tional box and random realization of the initial conditions. 
This stems for the fact that the LG constitutes a quasi-linear 
object, far from being in virial equilibrium. LG-like objects 
are delicately defined to match the observed LG and the 
small dynamical differences introduced by the A term are 
likely to prevent a full correspondence between objects in 
the different models. 

Some of the simulations used here are constrained ones, 
namely their large scale structure is constrained to repro- 
duce the observed cosmological neighborhood. No attempt 
is made here to select only LG-like objects that reside in the 
'correct' position within the cosmic web that constitutes the 
local neighborhood, as this would result in a very poor statis- 
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MODEL 


CONSTRAINED 


BOX [h-^Mpc] 




Qa 


h 


0-8 


N 


mDM[h-^MQ] 


ACDM 


yes 


64 


0.30 


0.70 


0.70 


0.90 


256=* 


1.3 X 10^ 


ACDMu 


no 


64 


0.30 


0.70 


0.70 


0.90 


256=* 


1.3 X 10^ 


ACDMhr 


yes 


64 


0.24 


0.76 


0.73 


0.75 


1024^ 


1.6 X lO'' 


ACDM160 


yes 


160 


0.24 


0.76 


0.73 


0.75 


1024^ 


2.5 X 10* 


OCDM 


yes 


64 


0.30 





0.70 


0.90 


256=^ 


1.3 X 10" 


OCDMu 


no 


64 


0.30 





0.70 


0.90 


256^ 


1.3 X 10" 


SCDM 


yes 


64 


1.0 





0.50 


0.7 


256=^ 


4.4 X 10" 



Table 1. Description of the set-up and cosmological parameters used for the different simulations: constrained or random, computational 
box size, matter density (Cm), cosmological constant (Oa), Hubble's constant (h), initial power spectrum normalization (eg), number 
of particles {N) and mass of the DM particle (m^jjif ) 



Criterion 


ACDM 


ACDMu 


ACDMhr 


ACDM160 


OCDM 


OCDMu 


SCDM 


Pairs 


13 


20 


12 


131 


11 


15 


24 


Singles 


52 


65 


43 


478 


38 


60 


78 



Table 2. Number of LG candidates found in each simulation following the Pairs and Singles criteria. 



tics. Rather, objects are defined as LG-like, Pairs or Singles, 
regardless of their location in the computational box, and 
the flow field around these objects is studied. 



5 ANALYSIS OF HUBBLE FLOW WITHIN 
THE LOCAL VOLUME 

The main purpose of the present section is to study the 
dispersion of the peculiar v elocities around LG- like objects 
selected in the manner of iMaccio etall (|2005l l. The flow 
around each object is studied using the RMS estimator of 
aH, so as t o allow a comparison with the analysis of the 
CNG data jTikhonov fc Karachentse^l2006l) . The analysis 
consists of three parts: 

(a) The examination of the Hubble diagram of the flow 
around LG-like objects. 

(b) The statistical distribution of an of the various models. 

(c) The dependence of an on the mean density of the LV 
around each object. 

In all cases we take a spherical shell of (1 — 7) Mpc as the LV 
of each object and, unless otherwise stated, an corresponds 
to all objects within the LV. Like with the actual data, the 
inclusion or omission of halos within 1.0 Mpc distance leaves 
the value of aH virtually unchanged. 

The six panels of Fig. [2] present the combined Hubble 
diagram of 10 randomly chosen LG-like objects of the differ- 
ent models. The solid line corresponds to the unperturbed 
Hubble flow and the upper and lower dashed lines corre- 
sponds to Hor ± an, where (jh is the median value of an 
for each simulation. The individual p oints with the error 
bars r epresent the value of (jh{R) from iTikhonov fc Klvpiril 
(|2008t) observational data. The results from ACDMhr sim- 
ulation are not shown in Figures [2] [3] and |4] since they are 
very similar to the other ACDM simulations. 

Table 13] presents the mean <jh of the full LV taken over 
the Pairs and the Singles in each of the simulations. The 
fractional cumulative distribution, riian) = A'^(< oh)/Nt, 



is presented in Fig. [31 where A*'(< an) is the number of 
LG-like objects colder than an and A'^r is the total number 
of LG objects. The cumulative distribution is shown for all 
models and for both the Pairs and the Singles. In each frame 
the vertical line indicates the observational value of an = 
90.4 km s~^ and the numbers give the fraction of LG-like 
objects with an smaller than the observed value of the Pairs 
and the Singles (in parentheses). 

Fig. |3]shows that there is a large scatter in the distribu- 
tion of (T£f and that for all the models considered here there 
is a non-negligible probability to find LG-like objects with 
(Jh equal or smaller the one corresponding to the actual LG. 
However, there is another dynamical variable that controls 
the dynamics around these objects and that is the mean den - 
sity within the LV. This was shown bv lMaccio et al.l (|2005l ) 
and it is clearly reproduced by the analysis of the models 
considered here. Fig. [J] shows the scatter plot of uh vs. the 
mean matter density within the LV normalized by the mean 
cosmological density. Both Pairs and Singles LG-like ob- 
jects are studies and the results are virtually the same in 
the range of Pdm/pdm ^ 2, where pdm and pDM are the 
mean DM density within the LV and the cosmological mean 
DM density. In the simulations the DM density is very easily 
measured but this is not a directly observable quantity. Ob- 
servationally the galaxy distribution within the LV is used 
to infer a local density. Short of having numerical simula- 
tions with full galaxy formation treatment the density of 
DM halos is used here as a proxy to the galaxy distribution. 
The lower panel of Fig. 3] presents the dependence of oh on 
the density of mass which belongs to halos within the LV 
normalized by the mean mass density of all the simulated 
volume. In both panels of Fig. [4] the horizontal line denotes 
the observational estimate of an. 



6 RELAXING THE LOCAL GROUP CRITERIA 

The LG appears to be a very typical small group of galax- 
ies, whose main dynamical characteristics are summarized 
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Figure 2. Combined Hubble diagrams of 10 randomly chosen LG-like objects in the different simulations. The solid lines represent the 
Hq r Hubble flow and the dashed ones show the median of ajj (see Table [3j. The individual data points (in red) correspond to the moan 
radial velocity and the ajj dispersion in the (1 — R) Mpc distance cut of the Tikhonov & Klypin (2008) data. The upper panel presents 
the LG-likc objects selected using the Pairs criterion and the lower ones correspond to those LG's selected with the Singles. 



Criterion 


ACDM 


ACDMu 


ACDMhr 


ACDM160 


OCDM 


OCDMu 


SCDM 


Pairs 
Singles 


102 ± 40 
70 ± 41 


91 ± 47 
80 ± 44 


90 ± 60 
69 ± 42 


114 ± 58 
87 ± 46 


80 ± 25 
60 ± 22 


81 ± 30 
60 ± 30 


117 ± 33 
74 ± 38 



Table 3. The median of au (in km s ^ units) and the standard deviation from the median for all the LG candidates in each simulation 
selected with Pairs and Singles criteria (sec Fig. [2ll. 
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Figure 3. Cumulative number of LG candidates, from Pairs (solid lines) and Singles (dashed lines) criteria, with a cth lower than a 
given value, for the different simulations. The vertical line shows the value of the observational cth = 90.4km s~^. The numbers on each 
frame are the fraction of LG-like objects with cr^f below this observational value for the Pairs and Singles (in parentheses). 





Distance between members (Mpc) 


No neighbours (Mpc) 


Virgo clusters 


(JH ( km s ^ ) 


a 


1.0 


3 


only one 


102 ± 40 


b 


singles 


3 


only one 


70 ± 41 


c 


1.5 


3 


only one 


85 ± 47 


d 


1.0 


no constrain 


only one 


189 ± 100 


e 


1.0 


5 


only one 


104 ± 14 


f 


1.0 


3 


no constrain 


113 ± 112 


g 


1.0 


3 


one or more 


150 ± 74 



Table 4. The median of cru and the standard deviation from the median for the candidates in the ACDM simulation for different 
selection criteria (see text). 



in ^ These are formulated in terms of the five criteria for 
the selection of LG-like objects. These dynamical properties 
are common in the Universe, yet the particular dynamical 
configuration of the LG determines an- It is interesting to 
see how these properties affect the flow field around the LG 
and for that purpose we have relaxed some of the criteria 
of fjH reselected ensembles of mock LGs and studied the cu- 
mulative histogram of an for the various selections (in the 
ACDM simulation). This is shown in Tableland Fig. O 
where the median and the cummulative fraction of LG-like 
objects of the constrained ACDM simulation, respectively, 
is shown for the following cases: 

(a) The Pairs, which are taken as a benchmark. 

(b) The case of Singles. 

(c) The Pairs criteria, but assuming a distance between LG 
members lower than 1.5 Mpc. 

(d) No constraint is imposed on the lack of nearby galactic 
size halos. 



(e) There must not be objects similar to LG members 
within a distance of 5 Mpc. 

(f) No constrained is imposed on the existence of a Virgo- 
like cluster. 

(g) There must be one or more Virgo-like clusters. 

The Pairs selection is taken here as a benchmark and it 
is represented here for the sake of completeness. The largest 
departure from the benchmark Pairs is obtained by remov- 
ing the constraint on nearby LG-like objects (case d). Table 
[4] and Fig. [5] clearly show the role of clusters in heating the 
flow. A nearby Virgo-like cluster exerts a tidal field in the 
vicinity of the LG-like object, resulting in a shear flow which 
contributes to the anisotropic component of the velocity field 
and thereby heating the flow. The benchmark case (a) has 
only one Virgo-like cluster while in case (g) there can be one 
or more clusters. As expected, the flow for objects found in 
case (g) is significantly hotter than for those found follow- 
ing case (a). The removal of the Virgo constraint increases 
somewhat the median but it more than doubles the stan- 
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Figure 4. A scatter plot of the relation between an and the mean matter density (normalized by the mean cosmological density) 
within the LV {R = 7Mpc) for Pairs (solid points) and Singles (crosses). Power law fits to the ajj - density scatter arc shown for the 
Pairs (solid line) and SingZes (dashed line). The horizontal line indicates the observational cr^f . In the upper panel the matter density is 
calculated from all DM particles inside the sphere. In the lower panel, the density is calculated from DM particles belonging to halos 
only. 



dard deviation (case f). In such a case mock LGs can either 
have rich clusters within their LV, hence with very high gh, 
or can reside in low density regions far away from clusters 
and consequently have very low gh- Similarly, allowing for 
nearby MW size halos heats the flow and leads to many hot 
LG-like systems (case d). The analysis shows that the most 
significant selection criterium is that of the nearby halos (d) . 
By omitting this criterion only roughly 15% of the LG-like 
objects have gh smaller than the observed value. 



7 MASS RESOLUTION 

The numerical simulations used here span a mass resolution 
that ranges over more than two orders of magnitude. The 
majority of the simulations are done on a very coarse grid 
of 256'^ and have low mass resolution. In fact, the LG-like 
objects of the low resolution ACDM and OCDM are made 
of about a thousand DM particles and the minimal mass of 
DM halos is ^ 3 X lO^°h"^M0. The faint nearby galaxies 
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Figure 5. Fractional cumulative distribution of for the LG 
candidates in the ACDM simulation found with the different cri- 
teria shown is Table |4] The black lines correspond to a (solid), 
6 (dotted) and c (dashed), the red ones to d (dotted) and e 
(dashed), and the blue ones to / (dotted) and g (dashed) cri- 
teria. The vertical line shows the value of the observational an- 



are presumably much less massive than this lower limit and 
the question arises as to what extent the present results are 
affected by the lack of resolution. 

To address the issue of mass resolution we have ana- 
lyzed the ACDMhr simulation and compared it with the 
ACDM one. This is a high resolution {N = 1024^) ACDM 
(WMAP3) constrained simulation in the same 64 hT^ Mc 
volume. This simulation differs from the low resolution 
ACDM one in terms of resolution, cosmological parameters 
and hence the power spectrum. It constitutes a different re- 
alization of the random field used to set the constrained 
initial conditions. An ensemble of LG-like objects has been 
constructed and oh has been calculated for all the DM ha- 
los in the corresponding LV (Figure |6] red curves) and for 
DM halos above the mass limit of the ACDM simulation 
(« 3 X 10^° h"^M0; blue curves). These are compared with 
the cumulative distribution of the low resolution ACDM sim- 
ulation. Figure [S] shows very clearly that, within the range 
studied here, the distribution of gh does not depend on the 
mass resolution. 



8 THE LOCAL VOLUME AND BLUESHIFTED 
GALAXIES 

A striking manifestation of the coldness of the local flow 
is the absence of blueshifted galaxies in the local neigh- 
borhood. In particular the CNG catalog shows only one 
blueshifted galaxy in the distance range of (1 — 7) Mpc, 
and it is flagged as a probable member of the Virgo, i.e. 
its distance might be erroneous. The absence of blueshifted 
galaxies provides only a qualitative measure for the coldness 



Figure 6. The fractional cumulative distribution (r)((T£f)), of LG 
candidates for simulations with different mass resolutions: The 
black curves correspond to the ACDM simulation and the red 
ones to the ACDMhr. The blue curves correspond to ACDMhr, 
but using only halos within the LG's in the same mass limit than 
those in the ACDM simulation. The solid lines correspond to 
the LG-Iike objects selecte by the Pairs and the dashed ones to 
those selected by the Singles criteria. The vertical line show the 
observational au = 90.4km s^^. 



of flow, yet it consti tutes a very selectiv e criterion for find- 
ing LG-like objects jSchlegel et al.|[l994l ') . The implications 
of such a selection on the coldness of the fiow around the 
selected groups are studied here. 

A word of cautious is due before applying the con- 
straint of no blueshifted 'galaxies' to the LG-like objects. 
Here, the local velocity field is traced by DM halos. It 
follows that the number of blueshifted objects might de- 
pend on the mass resolution of the simulations. We fol- 
low here iTikhonov fc KlvpinI { 20081 ) and assume that in 
the ACDM model DM halos with Vdrc > 35 km s~^ serve 
as proxy for observed luminous galaxies in the LV. Such 
halos are well resolved in the two high resolution simula- 
tions to which the no blueshifted 'galaxies' selection is ap- 
plied. Table [5] lists the results of an analysis applied to the 
ACDMhr and ACDM160 simulations. The analysis is ap- 
plied to DM halos with Vdrc > 35 km s^^. The analysis is 
applied to all the LG-like objects and then to groups that 
contain no blueshifted galaxies in the shell corresponding 
to (x - 7) Mpc, where a; = 2, 3 and 5 Mpc. The Table lists 
the number of LG-like objects and the median and standard 
deviation of (Jh of the selected objects. Figure [7] shows the 
cumulative histogram of an for the ACDMhr and ACDM160 
simulations. The figure presents the benchmark case of all 
pairs and the objects with no blueshifted Vdrc > 35 km s~^ 
DM halos in the 5 -7 Mpc and 3-7 Mpc sheUs. Figure [H] 
shows a scatter plot of an against the mean density within 
the LV, for the same objects considered in Figure [T] 

As expected, both Tableland Figure [7] show that the 
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constraint of no having blueshifted galaxies implies LG-like 
objects with cold environments. Both high-resolution sim- 
ulations have yielded objects with no blueshifted galaxies. 
As can be deduced from Table (5] the closer to the LG we 
look for blueshifted galaxies, the less number of LG's pass 
the criterion. Given the uncertainties involved in the asso- 
ciation of Vcirc 35 km s~^ halos with nearby faint galax- 
ies, we think that the number of LG-like objects with no 
blueshifted galaxies cannot be used as a reliable statistic of 
the local environment. Yet, it provides another strong obser- 
vational evidence for the coldness of the local flow. Figure 
[S] provides further support to our earlier findings about the 
relation between gh and the mean density of the LV (Fig- 
ure|4ll. LG-like objects with no blueshifted 'galaxies' in their 
LV are cold, gh ^ 100 km s^^, and they resides in LV with 
Pdm/pdm ^ 2. 



9 DISCUSSION 

A suite of CDM simulations have been used to study the 
nature of the velocity field around Local Group - like ob- 
jects. The objects are selected so as to reproduce the main 
dynamical properties of the LG, which include the range of 
masses of the two most massive DM halos in the group, the 
proximity to a Virgo-like cluster and the absence of nearby 
massive halos. The flow fleld has been analyzed by means of 
aH , which is defined by the RMS value of the residual from 
the pure Hubble flow, of all galaxies within the Local Vol- 
ume, deflned to be the volume enclosed within a radius of 
7Mpc around the LG. The results are compared with the an 
calculated for the CNG data. Two important conclusions are 
reached here. The first, and the more important one, is that 
there is no 'coldness of the local flow' problem, but rather it 
is the relation between an and the mean local matter den- 
sity that needs to be addressed in the context of the LV. The 
other is that within the canonical CDM cosmological mod- 
els an depends only on the mean ma tter density and not on 
the cosmological constant. Recentlv. iTikhonov et all (|2009l ) 
have found a similar cold fiow around LG-like objects in a 
model with Warm Dark Matter (WDM) of 1 keV particle 
masses which leads to a cut-off in the power spectrum above 
fepeak = S.Th-^Mpc. 

The main result that has been found is that for all mod- 
els considered here there is a non-negligible fraction of the 
LG-like objects with an close to the observed value. The 
two ACDM simulations in the box 64h~^Mpc have about 
(30 - 45)%, the ACDM160 has 25%, the two OCDM simu- 
lations have roughly (50 — 65)% and even the SCDM model 
has 17%, of their (Pairs) LG-like objects colder than the ob- 
served value. However the scatter in an is not random but 
it correlates strongly with the mean matter density within 
the LV, with higher an expected for higher density. In the 
canonical ACDM case a an ^ 100 km s^^ implies that the 
mean LV density is less than twice mean cosmological den- 
sity. This is also consistent with the requirement that there 
are no blueshifted galaxies in the LV around the LG. It 
follows that the old 'coldness of local fiow' problem reap- 
pears as a potentially new problem, namely the local den- 
sity - coldness relation. The rece nt compilation of the abun- 
dance of galaxies in the LV of iTikhonov fc KlvpinI (|2008l ) 
provides interesting results. These authors studied the lu- 



minosity function of the galaxies in the CNG data and es- 
timated the LV fractional density to be about 1.4 ± 0.17 
within 7? = 8Mpc. Using DM halos with Vdrc > 35 km s~^ 
as a proxy to luminous galaxies, the ACDM predicted den- 
sity enhancement in the LV stands in good agreement with 
the observed value. The realization that aH and the local 
density are closely connected provides another reminder to 
the fact that progress in understanding the local dynam- 
ics depends on a parallel advance in understanding galaxy 
formation. 

It has been speculated that the coldness of the local fiow 
is a manifestation of the dark energy, which dominates the 
mass-energy density of the universe |Barvshev et al. I l200ll . 
IChernin et 'aLll2004l . l2007l . lMacci6 et al.ll2005h . The basic ar- 
gumentation behind that claim is that very locally the grav- 
itational field is supposed to be dominated by the cosmolog- 
ical constant and hence the induced Hubble fiow is expected 
to be cold. This claim has been refuted by iHoffman et al.l 
(|2008l ) who showed that within a local volume defined by 
0.75 ^ i? ^ 2 Mpc the fiow is dominated by the DM, re- 
gardless of the existence of a cosmological constant. Namely, 
the OCDM and ACDM models yield virtually the same an 
within that volume. The current paper exte nds the numeri- 
cal ex periments and reinforces the results of iHoffman et all 
(|2008l ). That earlier work is extended to cover more simu- 
lations of higher resolution and the LV is defined to extend 
over 1 _R ^ 7 Mpc. 

The current work substantiate and giv es justification 
to the theoretical arguments suggested by IHoffman et al.l 
(|2008l ) as to why dark energy does not dictates the local 
dynamics. The arguments on the role of the dark energy 
are all based on approximating the local dynamics by the 
monopole term, i.e. a gravitational field induced by an iso- 
lated point-like object. This would have implied a very tight 
correlation between an and the mean overdensity within 
the LV, contrary to the large scatter manifested by Figure 
[T] The arguments advocating the cooling of the local fiow by 
the dark energy completely ignore the environmental effects 
on the flow. These effects are very clearly demonstrated in 
Table |4] which shows that aH nearly doubles by removing 
from the benchmark criteria the constraint on the lack of 
MW-like halos within a sphere of 3 Mpc. The proximity of 
Virgo-like halos also changes the value of an ■ It follows 
that the dynamics of quasi-linear objects such as the LG 
cannot be properly modeled without a detailed accounting 
of environmental effects. 

The present paper reac hes radically differen t conclu- 
sions than those expressed in lMaccio et al.l |200^ concern- 
ing the role of the cosmological constant in shaping the dy- 
namics of the LV. Yet, a close inspection shows that there is 
no technical disagreement between the two papers. We were 
careful to follow th e selection criteria of LG-like objects of 
iMaccio et aH (jiooi) and hence their ACDM simulation can 
be directly compared with the present ones. In spite of the 
somewhat different value of the physical parameters of the 
ACDM model used in the two papers the scatter of an and 
its dependence of th e local density are i n close agreement. 
A careful reading of lMaccio et al.l (|2005l ) reveals that their 
statement concerning the role of the A term is based on the 
comparison these authors did with the OCDM simulation of 
ICovernato et all (|l997h . 

A careful comparison with the an - Sp/p relation of 
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no negative velocity halos ACDMhr ACDM160 





LG candidates 


median 


LG candidates 


median 


all 


12 


97 ± 61 


131 


115 ± 58 


5-7 Mpc 


7 


85 ± 15 


46 


90 ± 31 


3-7 Mpc 


4 


85 ± 16 


15 


95 ± 31 


2-7 Mpc 


1 


69 ± - 


6 


99 ± 19 



Table 5. The case of LG-like objects with no blueshifted galaxies: The left column presents the selection criteria of the LG-like objects: 
all - the benchmark Pairs, the following raws corresponds to LG-like objects with no blueshifted DM halos with Vdre > 35 km s~^, 
taken as a proxy to observed galaxies, within the respective shell. Each raw provides the number of such groups and the median of cru 
and its standard deviation measured within their LV. The analysis was applied to the two high resolution simulations, ACDMhr and 
ACDM160. 
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Figure 7. The Fractional cumulative distribution of LG candidates (»?(o'£f )) for the two high resolution ACDMhr (left) and ACDM160 
(right) simulations, the solid line corresponds to the Pairs criterion, the other two lines represents groups with no blueshifted DM halos 
with Vcirc > 35km s~^ within 5-7 Mpc (red-dotted) and 3-7 Mpc (bluc-dashed) shells. The vertical line shows the observational 
(Jh = 90.4km s"^ 



the L G candidates in the OC DM simulation shown in Fig. 
11 of Governato et al. (|l997l ). and its presentation in Fig. 
2 of iMaccio et al. I I2005I . reveals that it can be reproduced 
in the present OCDM simulations by omitting the Virgo 
constraints on the selection of the LG-like objects. Given 
the present much better resolution, the n umber of our LG 
candid ates is considerable larger than in iGovernato et al.] 
j 19971 ). in particular in t he densest enviro nments. The lin- 
ear fit shown in Fig. 11 of lGovernato et al] is biased towards 
the LGs located in 5p/ p ^ 1, where most of their LGs are 
found. Restricting our fit to LGs with 5p/p ^ 1, a quite 
similar fit emerges. This indicates that the old and present 
OCDM simulations equally reproduce, to within their reso- 
lution limits, the same kind of LG candidates. To summa- 
rize, by neglecting the Virgo constraints on the selection of 
LG-like objects and by restricting the analysis to less dense 
objects the old OCDM results a re recovered. This might hint 
that the OCDM an analysis of lGovernato et al] (|l997l ) was 
performed with respect to their 'cat2' catalog rather than 
the stated 'cat3' catalog, in which the Virgo constraints were 
imposed. 

No substantial differences are found between the con- 
strained and unconstrained simulations. The constraints im- 
posed here a re affecting the stru cture on scales larger than 
« Sh-^Mpc (|Klvpin et al.ll2003l ). a scale that coincides with 
the LV. This implies that the behavior of gh predominantly 



depends on the internal dynamics within the LV and less on 
the one induced by larger scales. This also gives further sup- 
port to the claim that, apart from the particularities of the 
very nearby structure, the universe around us constitutes a 
very typical realization of the CDM cosmogonies. The par- 
ticular dynamical attributes of the LG are its mass range 
of about 10^^ h~^M0, its binary nature, the proximity to 
the Virgo cluster and the absence of similar nearby objects. 
These have been used as the selection criteria of LG-like 
objects in the simulations. In ^these criteria have been re- 
laxed, while keeping the mass to be roughly the LG mass. 
One expects that the proximity of the Virgo cluster heats the 
flow field in the vicinity of the LG, due to its tidal field which 
induces non-isotropic motions around the LG. Hence the re- 
laxation of the condition for a nearby Virgo-like cluster leads 
to colder LG-like objects. Similarly, in LG-like objects with 
a single MW-like halo the gravitational field is dominated 
by a monopole term and this leads to a colder fiow than in 
binary LG-like systems. On the other hand the largest, yet 
opposite, effect is introduced by allowing MW-like halos to 
be close to the LG. The presence of such nearby halos heats 
the flow considerably. It follows that the main ingredient 
that is responsible for the coldness of the local flow is the 
relative isolation, and therefore also the low mean local den- 
sity, of the LG. Furthermore, in the CDM cosmologies the 
coldness of the local flow depends more on the properties 
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Figure 8. Scatter plots of ajj vs. the mean matter density, normalized by the mean cosmological density, within the LV's found in the 
high resolution simulations, ACDMhr (left) and ACDM160 (right). All symbols correspond to LV's around LG-like objects found with 
the Pairs criterion in both simulations. The filled circles represent LV's with no blueshifted DM halos with V^irc > 35 km s~^ within 
the 5-7 Mpc shell (blue) and within 3-7 Mpc shell (red). The crosses are for the rest of the LV's. 



of the LG than on the possible presence of a dark energy 
component. 
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